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We report the first three-particle coincidence measurement in pseudorapidity (A77) between a high 
transverse momentum (p±) trigger particle and two lower p± associated particles within azimuth 
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|A0|<O.7 in = 200 GeV d+Au and Au+Au collisions. Charge ordering properties are ex- 

ploited to separate the jet-like component and the ridge (long-range Arj correlation). The results 
indicate that the particles from the ridge are uncorrelated in Arj not only with the trigger particle 
but also between themselves event-by-event. In addition, the production of the ridge appears to be 
uncorrelated to the presence of the narrow jet-like component. 

PACS numbers: 25.27.Gz 



Di-hadron coincidence measurements provide a pow- 
erful tool to study the properties of the medium cre- 
ated in ultra-relativistic heavy-ion collisions. The ob- 
servation of the long range pseudorapidity correlation in 
central Au+Au collisions lj, called the ridge @], where 
hadrons are correlated with a high transverse momentum 
trigger particle in azimuth (A^>~0) but extended to 
large relative pseudorapidity (A77), has generated great 
interest. Various theoretical models are proposed to ex- 
plain this phenomenon, including (i) longitudinal flow 
push H , (ii) broadening of quenched jets in turbulent 
color fields [4|, (iii) recombination between thermal and 
shower partons H, (iv) elastic collisions between hard 
and medium partons (momentum kick) and (v) par- 
ticle excess due to QCD bremsstrahlung or color flux 
tube fluctuations focused by transverse radial flow 0-QjJ . 
Models (i)-(iv) attribute the ridge to jet- medium interac- 
tions: particles from jet fragmentation in vacuum result 
in a peak at At7~0 and those affected by the medium 
are diffused broadly in Arj forming the ridge. Model (v) 
attributes the ridge to the medium itself, and its correla- 
tion with high p± particles is due to the transverse radial 
flow. 

Despite very different physics mechanisms, all mod- 
els [3l4l 1 1 give qualitatively similar distributions of cor- 



related hadrons with a high-pj^ trigger particle. Some 
of these model ambiguity can be lifted by 3-particle co- 
incidence measurements. We analyze the hadron pair 
densities from 3-particle coincidence measurements in 
(A?7i,A?72), the pseudorapidity differences between two 
associated particles and a trigger particle. We exploit 
charge combinations in an attempt to separate the jet- 
like and ridge components and study their distribu- 
tions, without assuming the A77 shape of the ridge. 
Jet fragmentation in vacuum should give a peak at 
(A?7i,Aj72)~(0,0), while particles from the ridge would 
produce structures that depend on its physics mecha- 
nism. Correlation between particles from jet fragmenta- 
tion and the ridge would generate horizontal or vertical 
stripes (At^i^O or A?72^0) in the 3-particle coincidence 
measurement. 

Results are reported for minimum bias d+Au, pe- 
ripheral 40-80% and central 0-12% Au+Au collisions at 
JWn = 200 GeV from the STAR experiment 12]. The 
40-80% data are from the minimum bias sample, and the 
0-12% data are triggered by the Zero Degree Calorime- 
ters (ZDC) in combination with the Central Trigger Bar- 
rel (CTB). This analysis uses 6.5 x 10 6 d+Au events 



taken in 2003, and 6.0 x 10 6 peripheral and 1.9 x 10 7 
central Au+Au events taken in 2004. The data are an- 
alyzed in finer centrality bins for Au+Au collisions [l3[ 
and are combined for better statistics. 

The reconstructed event vertex is restricted within 
|-Zvtx|<30 cm along the beam line from the center of the 
STAR Time Projection Chamber (TPC) [14], which sits 
in a uniform 0.5 T magnetic field. The data were taken 
with both magnetic field polarities. The trigger and as- 
sociated particles are restricted to |??|<1 and their p± 
ranges are 3<p^<10 GeV/c and l<p±<3 GeV/c, re- 
spectively. The correlated single and pair densities with 
trigger particle are corrected for the centrality-, p±-, <fi- 
dependent reconstruction efficiency for associated parti- 
cles and the (^-dependent efficiency for trigger particles, 
and are normalized per corrected trigger particle. 

Due to the high TPC occupancy of Au+Au events, 
track pairs close in rj and <f> can be merged and recon- 
structed as single tracks. This results in deficits in pair 
density at A?7^0 and at small, but non-zero, A<f) whose 
value depends on p± , charge combination and magnetic 
field polarity. To reduce this effect, we apply cuts to ex- 
clude close track pairs in real and mixed events. Losses 
due to those cuts are compensated for by the accep- 
tance correction obtained from mixed events. To en- 
sure the mixed events have similar characteristics as the 
real events, we mix events from the same centrality bin 
without requiring a trigger particle and with the same 
magnetic field polarity and nearly identical z v tx position, 
refered to hereon as inclusive events. 



Au+Au 0-12%, 



(b) Au+Au 0-12%, |A(M<0.7 




1 2 3 

Ai> (radians) 



FIG. 1: Correlated hadron distribution in (a) A<j> (|»j|<l), 
and (b) Arj (\A(f>\<0.7) with a high-px trigger particle in 0- 
12% Au+Au collisions for 3<p ( L t) <10 GeV/c and Kp^O 
GeV/c. The ZYAl-normalized flow background is shown in 
(a) by the curve. The Arj distributions in (b) are background 
subtracted and corrected for Arj acceptance, and are for like- 
and unlike-sign pairs separately. The curves in (b) are Gaus- 
sian fits. Errors are statistical. 



4 



Figure QJa) shows the hadron Acf> distributions rel- 
ative to the trigger particle in 0-12% Au+Au col- 
lisions. Also shown is the background B(A<fi) — 
aF{A4>) f„ Bi nc (Ari, A(f>)dArj where i?i nc is constructed 
by mixing a trigger particle with associated particles from 
a different and inclusive event. The flow contribution 

F(A<p) = 1 + 2i4*M a) cos (2A</>) + 2v [ l ] vf cos (4A</>) 

(1) 

is added to mixed events using the measured, rj- 

independent, v 2 [IH and a parameterization of 

«4=1.15t;| |13|. A normalization factor, a, is ap- 
plied to match the distribution in O.8<A0<1.2, as- 
suming zero yield at A<t>^l radian (ZYA1) [l]. The 
near-side (|A0|<O.7) correlated hadron yield in A-q is 
Y(Ar)) = Y(Ari) - B(Arj), where Y(Arj) and B(Arj) = 
a J° Q 7 7 -Binc(Af?, A(j))F(A(j))dA(j) are the signal and back- 
ground distributions, respectively. Figure[ljb) shows the 
Y(At]) distribution, after 2-particle Arj acceptance cor- 
rection, for the like- and unlike-sign trigger-correlated 
particle pairs. Jet-like peaks at Arj^O are observed, 
atop a broad, charge- independent pedestal (the ridge). 
A Gaussian fit to the peak yields o=0.50±0.04 for the 
like-sign and <t=0.41±0.01 for the unlike-sign pairs in 
Arj. 

All triplets of one trigger particle and two associated 
particles from the same event within |A0i j2 |<O.7 are 
analyzed. Combinatorial background B\ (or B2) arises 
where only one (or neither) of the two associated particles 
is correlated with the trigger particle besides flow corre- 
lation [l6| . The former cannot be readily obtained from 
the product of the event averaged Y(Avi) and B(Arj), be- 
cause of the varying Ar) acceptance from event to event. 
Instead, we construct B\ by mixing trigger-associated 
pairs from the real event with a particle from a different 
and inclusive event, namely, 



Bi = 



aY(A m )B inc (A m )(F^ 2 \Afc) + 



F^^A^-A^ + F'-l) + (lo2) 

2 a 2 B i „c(Ar, 1 )B i „c(A? 72 )(F( t ' 1 )(A0 1 )-|- 
F^(A4> 2 ) + F^fA^ - Afa) + F' - 2) 



(2) 



Here the last term is constructed by mixing the trigger 
particle with two different inclusive events to remove the 
uncorrelatcd part in the first two terms, and 



F' = 2v% ) v'£ ) vf ] cos(2A0i - 4A</> 2 ) 
+ 2v < i ) vf ) v i l ) cos(4A0i - 2A</> 2 ) 
cos(2A0i + 2A</> 2 ). 



2v. 



(1)„,(2) (t) 



(3) 



The flow terms [16] in (. . .) are added in because they are 
lost in the event-mixing; their averages are taken within 
|A0i j2 |<O.7. The superscripts represent the u 2 and v\ for 



trigger and associated particles. To increase statistics, 
we mix each trigger particle with ten different inclusive 
events. 

The second background (£> 2 ) is constructed by mixing 
a trigger particle with associated particle pairs from in- 
clusive events thereby preserving all correlations between 
the two associated particles (denoted by £g>) [l6j |: 



B 2 = a 2 6^B inc (A?7i) ® B inc (Ai] 2 
(F^A^+F^fA^ + F'-l). 



(4) 



The factor a 2 b scales the number of associated hadron 
pairs in the inclusive event to that in the back- 
ground underlying the triggered event: b = ((N(N — 

where N denotes 



l))/<A0 2 Wd/((A(A - 1))/(A0 2 ) 
the associated hadron multiplicity [16| . If the associated 
hadron multiplicity distributions in both the inclusive 
event and the background are Poissonian, or deviate from 
it equally, then 6=1. We obtain b as follows. We scale 
the correlated hadron A77 distribution such that there 
would be no ridge in 1.0<| A?y|<1.8, and this gives a new 
value for a. We repeat our analysis with this new a, and 
obtain b by requiring the average correlated hadron pair 
density in 1.0<|A?7i. 2 |<1.8 be zero. We use the obtained 
b with the default ZYAl a to obtain the final 3-particle 
coincidence signal. The assumption in this procedure is: 



(N(N-l))/(Ny 



(N(N - 



1)>/(A0 2 

bkgd L J bkgd+ridgc 

(5) 

and is reasonable gauged from multiplicity distributions 
of inclusive and triggered events. The background- 
subtracted correlated pair density is corrected for 
3-particle Ayy-Af? acceptance, which is obtained from 
event-mixing of a trigger particle with associated par- 
ticles from two different inclusive events. We use ten 
pairs of inclusive events for each trigger particle in the 
mixing. 

The main sources of systematic uncertainty in our 
results are those in a, b and u 2 . These uncertainties are 
mostly correlated, therefore having insignificant effect on 
the shapes of our correlated density distributions. The a 
and b values for 0-12% Au+Au collisions are 0.998±£;[$ 2 
(syst.) and 0.99986tg;^^ 2 (syst.), respectively. The 
uncertainty on a is estimated by using the normalization 
ranges of O.9<A0<1.I and 0.7<A</><!.3. That on 
b is estimated by using the normalization ranges of 
1.8<|A?7|<1.2 and 1.2< | A77I <0.6. We note that the 
ridge is defined under the assumption of ZYAl in A0, 
by the factor a. Deviations of a from this assumption 
are not included in our systematic uncertainties. Such 
deviations (e.g. 3-particle ZYAM [l^]) do not introduce 
significant change to the shape of the ridge. 

The V2 systematic range used in our analysis is given 
by those from the modified reaction plane and 4-particle 
cumulant methods [ij and their average is used as 
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FIG. 2: Background-subtracted charge-independent (AAT) correlated hadron pair density in (a) minimum bias d+Au, (b) 
40-80% Au+Au, and (c) 0-12% Au+Au collisions for 3<p^<10 GeV/c and Kp^'o GeV/c. The results are for near-side 
correlated hadrons within A01.2 1 <0.7, and corrected for the 3-particle A77-A77 acceptance. Statistical errors at (At7i,A7/2)~(0,0) 
are approximately 0.033, 0.058', 0.084 for d+Au, 40-80% and 0-12% Au+Au, respectively. 



our nominal i> 2 . An additional systematic uncertainty 
arises from possible correlation of the ridge with the 
reaction plane which is not included in Eq. @. The 
estimated uncertainty from this source and that from 
i>2 are added in quadrature and referred to generally as 
flow uncertainty. 

Figure [JJa-c) shows the background subtracted charge 
independent (referred to as AAT) correlated hadron pair 
density (P) for minimum bias d+Au, 40-80% and 0-12% 
Au+Au collisions, respectively. The d+Au and 40-80% 
Au+Au results show a peak at (A771, A^ 2 )~(0,0), 
consistent with jet fragmentation in vacuum. A similar 
peak is also observed in 0-12% Au+Au collisions, but it 
is atop an overall pedestal. This pedestal is composed of 
the ridge particle pairs, and does not seem to have other 
structures in (A?7i,A?72)- To see this quantitatively, 
Fig. Eta) shows the average (P) for AAT as a function 
of R=^Ar)l + ~Kr%. The average density is peaked 
at R^0 and decreases with R for all systems. For 
d+Au and 40-80% Au+Au collisions the average density 
at R>1 is consistent with zero, indicating no ridge 
contribution. On the other hand, in 0-12% Au+Au 
collisions, the average denstiy drops more slowly and 
becomes approximately constant above R>1, indicating 
the presence of the ridge. 

Jet fragmentation has a charge ordering property, 
as shown at |A?7|~0 in Fig. W[h). The probability to 
fragment into three same-sign hadrons at our energy 
scale is negligible 
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same-sign hadron triplets may therefore be interpreted 
as ridge correlation. Thus, we analyze our data with 
same-sign triplets (A ± A ± T ± ) only, as well as with a 
same-sign associated pair and an opposite-sign trigger 
particle (A ± A ± T T ). The results are shown in Fig.Etb). 
Indeed, no jet-like component is apparent in A^A^T^. 
The A ± A ± T T result contains both jet-like and ridge 
components. The contribution from other charge com- 
binations, namely A A Z ^T ± , are simply the difference 
between AAT in Fig. EK) and {A ± A ± T ± + A ± A ± T^) 



o d+Au 

□ Au+Au 40-80% 
A-Au+Au 0-1 2% 




0.1 
0.08 
0.06 
0.04 
0.02 




(b) Au+Au 0-12% 

A A ± A ± T + 
▲ A^T* 

rA- 



0.5 



1.5 



R=\/(A^+A^ 



R=\/(AT^+Ati; 



FIG. 3: The average correlated hadron pair density per trig- 
ger particle as a function of R (a) for all charges in mini- 
mum bias d+Au, 40-80% Au+Au and 0-12% Au+Au colli- 
sions, and (b) for same-sign associated particles (A A ± T ± 
and A ± A ± T 11 ) in 0-12% Au+Au collisions. Systematic un- 
certainties are shown in the shaded boxes due to background 
normalization and in the solid curves due to flow. 



in Fig. E^b). We found this to be equal to twice the 
A ± A ± T^ contribution within errors. 

The ridge is very similar for like- and unlike-sign 
trigger-associated pairs at |A?7|>0.7 as shown in 
Fig. QTb), thus we expect the ridge contributions in 
the correlated pair density to be the same in all charge 
combinations. We verified this for large A77 correlated 
pair densities within our current statistics, as can be 
seen from Fig. Etb). Therefore the total ridge parti- 
cle pair density (P rr ) can be obtained as four times 
A ± A ± T ± . The remaining jet-like signal, the sum of 
jet-like correlated particle pairs (Pjj) and cross pairs of a 
jet-like and a ridge particle (Pj r ), can then be obtained 
by subtracting the total ridge from AAT. 

Figure BJa) shows the R dependence of the average 
(P rr ) and (Pjj)+(Pjr) in 0-12% Au+Au collisions. 
The ridge pair density is consistent with a constant 
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0.14±0.02 (x 2 /ndf=5.8/7). Gaussian fits indicate a best 
fit value er=2.1 (x 2 /ndf=4.8/6, solid curve) and <7>1.4 
(dashed curve) with 84% confidence level. On the other 
hand, the jet-like component is narrow with a Gaussian 
a=0.34toll (x 2 / ndf =0-8/6, dominated by statistical 
errors), comparing well to those from the correlated 
single hadron density. 

In order to investigate possible structures in the 




FIG. 4: The average correlated hadron pair density per trig- 
ger particle in 0-12% Au+Au collisions (a) for the jet-like and 
ridge components as a function of R, and (b) for the ridge as 
a function of £ within 7?<1.4. The solid curves are Gaussian 
fits. The dashed curve is a Gaussian fit with a fixed a—lA 
(see text) to the ridge data. The systematic uncertainities on 
the ridge data are shown in shaded boxes due to background 
normalization and in open boxes due to flow. 

ridge, we show in Fig. HJb) the average ridge particle 
pair density as a function of £=arctan(Ary2/A?7i) within 
i?<1.4. The data are consistent with a uniform distribu- 
tion in £ (x 2 /ndf=1.7/7). This suggests that the ridge 
particles are uncorrelated in Arj not only with the trigger 
particle but also between themselves. In other words, 
the ridge appears to be uniform in Arj event-by-event. 

Correlation between the jet-like correlated hadrons 
and the ridge would yield horizontal and vertical stripes 
in the correlated pair density in Fig. [21 resulting in a 
finite Pj r , a non-zero signal in Fig 0Ja) at large R. We 
found (P ir )=-0.001±0.030, averaged over the |A?7i|<0.7 
and |A?72|>0.7 region and its mirror region. On the 
other hand for the correlated jet-jet and ridge-ridge pairs, 



p jj ) ( p rr) = yj (0.081 ± 0.034) x (0.114 ±0.039)=0.096± 
0.026 where the averages are taken with |At7i )2 |<0.7 
and I A?7i 5 2| >0.7, respectively. The comparision between 
these two pair density magnitudes (whose systematic 
uncertainities are strongly correlated) suggests that 
production of the ridge and production of the jet-like 
particles may be uncorrelated. 

Our data qualitatively distinguish between some of 
the ridge models, (i) Longitudinal flow [1] would push 
correlated particles in one direction yielding a diagonal 
excess in A77-A77, disfavored by the present data, (ii) 
Turbulent color fields [H would generate a broad ridge in 
A77, which may however still be too narrow to reconcile 
with the width of our ridge pair density distribution. 



(iii) Recombination between thermal and shower par- 
tons Q should produce horizontal and vertical stripes in 
correlated pair density distribution which is disfavored 
by the data, and it does not have a mechanism for 
long range A77 correlations, (iv) The momentum kick 
model incorporates a broad ridge as input, but it should 
produce a much larger ridge on the away-side than 
on the near-side which is not supported by data fl9j |. 
and also may not describe other data such as the 
reaction plane dependence of the ridge in di-hadron 
correlations [13 ■ (v) QCD bremsstrahlung 0, H| or color 



flux tube fluctuations |9Ull| would yield a structure-less 
pair density [lot ] for the ridge as observed in our data, 
however the correlations between jet-like particles and 
ridge, as expected from these models, are not observed 
with our present sensitivity. Clearly more quantitative 
model calculations are needed to compare to the data 
reported here and elsewhere P, 0, Ho] to further our 
understanding of the ridge. 

In summary, we have presented the first 3-particle 
coincidence measurement in A?7-A?7 in minimum bias 
<i+Au, 40-80% and 0-12% Au+Au collisions at ^s^T/v = 
200 GeV. The p± ranges are 3<p ( | t) <10 GeV/c for the 
trigger particle and l<p^<3 GeV/c for both associated 
particles. A correlated hadron pair density peak at 
(A7yi,A772)~(0,0), characteristic of jet fragmentation, is 
observed in all systems. This peak sits atop a broad 
pedestal in 0-12% Au+Au collisions, which is composed 
of particle pairs from the ridge. We have exploited the 
charge ordering properties to separate the jet-like and 
ridge components. We found that same-sign associated 
pairs correlated with a same-sign trigger particle are 
dominated by the ridge. While the jet-like particle 
pair density is narrowly confined, the ridge is broadly 
distributed and is approximately uniform in Ary. A 
Gaussian fit in R to the average correlated pair density 
of the ridge yields ct>1.4 with 84% confidence level. 
Except for the correlations at Atfi^O, the particles from 
the ridge appear to be uncorrelated in Arj not only with 
the trigger particle, but also between themselves; they 
are uniform in our measured A 77 range event-by-event. 
No correlation is found between production of the ridge 
and production of the jet-like particles, suggesting the 
ridge may be formed from the bulk medium itself. 
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